Implications of a viscosity bound on black hole accretion 
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Motivated by the viscosity bound in gauge/gravity duality, we consider the ratio of shear viscosity (7/) to 
entropy density (s) in black hole accretion flows. We use both an ideal gas equation of state and the QCD 
equation of state obtained from lattice for the fluid accreting onto a Kerr black hole. The QCD equation of state 
is considered since the temperature of accreting matter is expected to approach 10 12 K in certain hot flows. We 
find that in both the cases 77/s is small only for primordial black holes and several orders of magnitude larger 
than any known fluid for stellar and supermassive black holes. We show that a lower bound on the mass of 
primordial black holes leads to a lower bound on 77/s and vice versa. Finally we speculate that the Shakura- 
Sunyaev viscosity parameter should decrease with increasing density and/or temperatures. 
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Accretion flows around black holes and neutron stars are 
very hot. In the geometrically thin (and optically thick, 
i.e., strong interaction between dense matter and radiation) 
regime, the temperature (T) of the flow could be ~ 10 7 K 
lM- The flows in certain temporal classes of micro-quasar 
GRS 1915+105 |0] are of this type. However, in the geo- 
metrically thick (and optically thin, i.e., weak interaction be- 
tween tenuous matter and radiation) hotter regime, we expect 
a two-temperature flow with the ion temperature (Tj) as high 
as 7 x 10 11 K and the electron temperature being much lower. 
The supermassive black hole system at the centre of our galaxy 
Sgr A* |3], (most of the temporal classes of) the X-ray binary 
Cyg X-l |Ql are of this type. The latter type of flows is radia- 
tively inefficient. 

In the existing literature, accretion flows are modelled us- 
ing an ideal gas equation of state (EoS) along with radiation. 
For hot optically thin flows, the underlying radiation mecha- 
nisms could be synchrotron, inverse-Comptonization of soft 
photons etc. The argument behind the choice of ideal gas is 
that the interaction energy between the constituent particles is 
small compared to their thermal energy. While this could be a 
viable option for cooler optically thick flows, the temperature 
of hot optically thin flows near the horizon of an astrophysi- 
cal black hole is of the order of the crossover temperature, T c 
(a temperature at low baryon density when the hadrons melt 
into a Quark Gluon Plasma (QGP)). At low densities (zero 
chemical potential), lattice QCD calculations |5J, |6[] suggest a 
crossover from the hadronic state to the QGP state at around 
1.5 — 2 x 10 12 K. At such high temperatures, there is expected 
to be copious pion production in a strongly interacting system. 
This QGP state is thought to have existed a few microseconds 
after the big bang and is being extensively studied at the Rela- 
tivistic Heavy Ion Collider (RHIC) and the Large Hadron Col- 
lider (LHC). Hence, the choice of ideal gas EoS seems sus- 
picious for the optically thin flows. Furthermore, ~ M & pri- 
mordial black holes [7] formed at T ~ 10 13 K (quark-hadron 
phase transition era) would indeed accrete QGP matter at high 
density ~ 10 18 gm/cm 3 . So the natural question arises: What 
EoS do we use for matter at such high temperature and/or high 
baryon density? At high temperatures and low baryon density, 
the answer is provided by lattice QCD. At high temperatures 



and high baryon density, the lattice methods have proved less 
useful and one would need to appeal to intuition drawn from 
models. 

Hydrodynamic simulations at the RHIC and LHC suggest a 
small value for the ratio of shear viscosity (77) to entropy den- 
sity (s) for the QGP state. Such a small value points towards 
a strongly interacting matter. String theory arguments and the 
famous gauge/gravity duality in the form of the AdS/CFT cor- 
respondence suggest a lower bound 18[] for r)/s: 
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While originally /i was thought to be unity JS], we know now 
that /i < 1 [9] and quite likely non-zero [10] (for recent re- 
views, see llllll ). For example, in the gravity model lfl2ll . 
/1 w 0.414. For the purpose of this work we assume that there 
is a bound and that /j, ~ 1. One of the main questions that we 
wish to answer is: How small is rj/s in black hole accretion? 
We can anticipate the answer to this question before getting 
into details. First, r\ is a combination of a molecular/physical 
viscosity and a turbulent viscosity. If the molecular viscosity 
dominates then for temperatures around or above T c , we ex- 
pect rj/s to be close to O.lh/ks- More precisely, we need to 
use the QCD EoS. If the turbulent viscosity dominates, then 
the ratio is much higher. However, as we will argue, close to 
the horizon where gravity effects on the flow are strong, or in 
the early universe when density was high, we expect a small 
ratio. 

There is a huge literature discussing accretion of optically 
thin hot matter around a Kerr black hole in full general rela- 
tivistic as well as pseudo-Newtonian frameworks (e.g. lfT3l - 
[l5l [nil ). For the present purpose, we will consider the 
one temperature advection-dominated accretion flow (ADAF) 
model around a Kerr black hole lfl4l [l5ll . The units are 
G = M = c = 1, with M being the mass of the black hole. 
In what follows we will work in terms of normalized tempera- 
tures where the normalization is done through m p c 2 /fcs with 
■trip being the mass of a proton. 

The equations of continuity and the stress energy conserva- 
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Here r is the Boyer-Lindquist radius, p the rest-mass density, 
He = H/r with H the half thickness of the flow, d the ra- 
dial velocity measured in the co-rotating frame, A the spe- 
cific angular momentum of the flow, V = 1 — 2/r + a 2 /r 2 
with a being the rotation parameter and M the rest-mass ac- 
cretion rate. The total energy density e = p + pTg(T) with 
g[T) = (4/(70 -l) + 15T)/(4 + 5T) and k = (e +p)/pis 
the relativistic enthalpy. The function / r is given by 
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with ^ = 1 + a 2 /r 2 + 2a 2 /r 3 , 7 2 = 1 + A 2 (l - d 2 )/{r 2 A) 
and fi± = ±(r 3 / 2 ± a)' 1 with 17 = 2a/(„4r 3 ) + A(l - 
■d 2 ) 1 / 2 ?) 1 / 2 / {^r 2 A 3/2 ) being the angular velocity. When 
complemented with the gas EoS p = pT, the equations © 
enable us to solve for T, p, A as functions of r. In (|2), j 
is a constant and should be treated as an eigenvalue for tran- 
sonic flows. The form of the equations enables us to rescale 
p — > Mp to set M = 1 in the numerics. To reinstate phys- 
ical units, for radial velocity we use -de, for specific angular 
momentum XGM/c, for density pMc 3 /G 2 M 2 and for tem- 
perature TmpC 2 /ks- The viscous stress tensor components 
entering the equations are WT and given by 
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with u r = -dV 1 / 2 j (1 — i? 2 ) 1 / 2 . The expression for a is very 
lengthy and can be found in the appendix of [ 14]. Parametriz- 
ing the turbulent viscosity based on the famous Shakura- 
Sunyaev prescription I 111 a nd accounting for the relativistic en- 
thalpy factor as in 014, 11511 we have 



r\ = anpCsHgr . 
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Here a is a constant typically taking values 0.01 — 0.1. The 
contribution from the physical (molecular) viscosity is small 
and the viscosity is turbulence induced. In the present work, 
we take a = 0.01 and 70 = 1.4444. The quantity r r is a 
relaxation time and is given by t,. = ri/(npc 2 s ). At this stage it 
is a phenomenological input. The presence of the relativistic 
enthalpy makes 77 proportional to e + p. Using the fact that the 
entropy density s = (e + p) jT « p(c 2 + c 2 )/T, we have 
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Following previous work ffl [l5ll . we reproduce the results 
using a shooting method shown in FIG.[T] We estimate a rough 
order of magnitude for 77/s for a = 0.999, when the minimum 
value of Hg is of the order of 0.3 near the horizon where k ~ 
0(1), c s - 0.4 and T - 0.3 for a = 0.01. This gives us 
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so that for stellar mass black holes, this ratio is of the order 
of 10 16 in units of h/kg, while for a supermassive black hole 
this is > 10 21 compared to the QGP value of 0.1! Only for 
primordial black holes with M/M Q ~ 10~ 16 will this ratio be 
comparable to the QGP value. In fact in order to be consistent 
with the viscosity bound, it appears that 
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Now, note the following curious fact. Had we been using a 
larger a ~ 1 we would have got 77/s ~ 10 17 M /MQU/ks- 
If we use the fact that in the present universe the surviving 
primordial black holes must have M > 10 15 gm we obtain 
77/s > Q.lh/k B \ 




FIG. 1: Variations of (a) Mach number, (b) temperature, (c) 77/s in 
units of h/kB, as functions of flow radius, when a = 0.999. 

FIG.[T]shows that close to the black hole event horizon T = 
0.25, corresponding to ~ 2.5 x 10 12 if, which is above T c . For 
QGP matter at this temperature we expect a small 77/s. Hence, 
the remarkably high value of 77/s for astrophysical black holes 
is quite puzzling. Even away from the event horizon, while T 
decreases (only an order of magnitude) below the crossover 
temperature, it is large enough to question the very large 77/s 
we find. One possible resolution to this puzzle is that at these 
temperatures the EoS being used is problematic and it is more 
appropriate to use the QCD EoS J5l|6[]. 
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The QCD EoS from lattice calculations is specified by the 
energy density e and pressure p. It is useful to express the 
result in terms of I(T) — e — 3p which is proportional to 
the trace anomaly. If we had conformal matter, this quantity 
would be zero. At high temperatures (T > 1.97 x 10 12 K), a 
useful way of fitting the lattice data is given by ^ 



I(T) 

rp 4 



J^2 



where d 2 = 0.2405GeI/ 2 , d 4 = 0.01355GeV 4 , ci = 
-0.0003237Gey 5 ,c 2 = 1.439 x 10~ 1A GeV w . A somewhat 
different fitting function is given in [5]. The crossover temper- 
ature in [5] works out to be around 150 MeV in comparison 
to what is used in 01, which is closer to 190 MeV. At low 
temperatures (< 1.97 x 10 12 K), 



= aiT + a 2 T 3 + a 3 T 4 + a 4 T 10 . 
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where ai = 4.654GeV r " 1 



a 2 = -879Gey- 3 ,a 3 = 
8081GeF- 4 ,a 4 = -7039000GeV r - 10 . The low temperature 
approximation is found using a matching procedure with the 
Hadron Resonance Gas model. From these we find 
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shown in FIG. 2. In order to go from lMeV to the normalized 
units we have to multiply by « 0.0011. 
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FIG. 2: Variations of (a) e/T 4 , (b)p/T 4 , as functions of temperature. 

In order to explain the QGP observables for heavy ion flows 
at RHIC/LHC, it has become apparent that a non-zero shear 
viscosity, -q, is needed. One way of characterizing r\ which 
is a dimensionful number, is to take its ratio with the entropy 
density. It turns out that the RHIC/LHC plasma has ill 611 
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Although 77 itself is large (in cgs units — 10 12 gm/cm s), its 
ratio with s is small. To put things into perspective, this ratio 
for water is around 2 orders of magnitude greater. Since for 
optically thin accretion flows (e.g. ADAF Hal . GAAF IU7I1 ) 
the temperature near the horizon can be of the order of 10 12 if , 
it is an interesting question to ask if rj/s becomes as small as 
in equation (fTTT >- 



The QCD EoS (see FIG. 2) cannot be valid everywhere 
since the low temperature behaviours of QCD matter and or- 
dinary matter are quite different, and far away from the black 
hole when flow is cooler, the accreting matter will presumably 
follow the ideal gas EoS. Hence, an accurate analysis would 
involve a knowledge of EoS that is valid everywhere, but this 
appears to be a very hard problem. However, we know that 
close to the horizon 77, (r < lOr^) flow must be transonic. 
Therefore, we start our computations for the QCD EoS with 
temperatures at an outer boundary that are similar to what we 
found using the ideal gas EoS. 

Before we start discussing the solutions, note that equations 
(|2]i-(8) exhibit the following scaling symmetry: (a) Rescaling 
p to p — > Mp and j — > j/M leaves the equations invariant, 
(b) Rescaling e, p and j to e — > (3e, p — >• f3p and j — > (3j leaves 
the equations invariant (f3 being any constant). Thus we do not 
need to carry the actual units of e,p, j. The relaxation time r r 
is taken to be r r — bri/(sT) [19]. Here b parametrizes our 
ignorance of the QCD coupling constant, influence of gravity 
etc. We choose b — 10 4 . The choice of b does not alter our 
conclusions. 

FIG. |3] shows the variations of Mach number and temper- 
ature as functions of radius. The variations ll20ll of r]/s for 
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FIG. 3: Variations of (a) Mach number, (b) temperature, as func- 
tions of flow radius. The solid dashed and dotted lines are for 
a — 0.999, 0.5 and respectively. 
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FIG. 4: Variations of rj/s in units of h/ks as a function of flow radius 
for the cases shown in FIG. [3] 



various values of a are shown in FIG.|4] This seems to suggest 
the same inequality as in equation @. 

So what have we learnt? In using the QCD EoS we have 
found the following: 
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• The temperature close to never reaches as high as 
what the ideal gas EoS dictates. 

• The value of rj/s is small only for primordial black holes 
and gives a result that is quite similar to what was found 
using the ideal gas EoS, for all black holes. 

• Using the fact that the mass of primordial black holes 
in the present universe must be greater than 10 15 gm, 
we find that rj/s > Q.lah/ks- Although for a < 1, 
r)/s < j^h/ kg, this is perfectly consistent with the fact 
that the original KSS bound can in fact be violated in 
holography |3. 

Let us conclude with some speculations. In order to apply 
the viscosity bound for black holes of mass O.IMq in the early 
universe around T c (which is in the radiation dominated era), 
one must consider an optically thick flow model lfll l2~lll . unlike 
the one we considered previously. In this quark-hadron phase 
transition era, the black holes must be accreting QGP matter. 
However, here we face the following conundrum. First, for 
QGP matter, the Reynolds number (Re) is 0(1). Hence, we 
do not expect turbulent viscosity to dominate in this era and 
we expect a small -q/ s. However, if we take a ~ 0.01 then 
rj/s comes out to be ~ O(l0 9 )H/kB in the inner accreting 
region in our computations. Second, the temperature in the 
inner region in this model |l|,|2l|] for a ~ 0.01 — 0.1 is only 
around 10 7 - 10 8 K. To get T > 10 n K for M = £M E dd JH, 
the value of a should be very small (< 10~ 14 ) which in turn 
leads to a small rj/s ~ £, 2 h/kB with Re ~ 10. Such a low a 
also leads to a large density which is consistent with the fact 
that the universe was very dense in this era. This seems to 
suggest that a should decrease with increasing density. Now, 
does it make sense to have rj/s > O(10 15 )ft/fcs near T c , 
which appears to be the case for stellar and supermassive black 
holes discussed above? Near the black hole, we expect grav- 



ity to make the flow sub-Keplerian and quasi-spherical, which 
is more "ordered", diminishing turbulence. This is supported 
by the fact that Re for QGP and the optically thin disk flows 
are comparable and as such if matter was in the form of QGP, 
turbulent viscosity would not dominate. A more natural so- 
lution is to demand that a becomes small near T c . Thus our 
conjecture is the following: 

The Shakura-Sunyaev viscosity parameter a should be a 
function of temperature and density such that it decreases with 
increasing temperature and/or density. However, low a would 
only support sub-Keplerian optically thin flows. Therefore, the 
inner accretion is necessarily sub-Keplerian. 

Finally let us point out a somewhat worrying aspect to hav- 
ing a large rj/s. First, we note that for the hydrodynamic ap- 
proach to be valid, it had better be true that the stress tensor is 
treated as an expansion in gradients of the velocity. In partic- 
ular = pg^ v + pnu^u 11 + , where is proportional 
to r) and involves gradients of the velocity. We have to ensure 
that the first order terms, t^ u , are much smaller than the ze- 
roth order terms, pg^ 1 pnu^u" . When we compute the ratio 
of these two terms along the flow, it turns out that the contribu- 
tion from t^ v is around 20 — 25% of that of the leading term. 
Hence, the numerics appear to be suspect since the subleading 
terms in the gradient expansion will be important. It may be 
worth investigating effects such as cavitation [23] in this con- 
text. It seems that lowering the viscosity parameter along the 
lines we have hinted at above may ameliorate this problem. 
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